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ABSTRACT 

We  treat  in  detail  positron  production  from  the  decay  of  radioactive  nuclei  produced  in  nuclear  reactions  of 
accelerated  ^He.  Because  of  their  large  cross  sections  and  low  threshold  energies,  these  reactions  can  significantly 
contribute  to  positron  production  in  solar  flares  with  accelerated-particle  compositions  enriched  in  ^He.  The 
addition  of  these  ^He  reactions  extends  earlier  calculations  of  positron  production  by  accelerated  protons  and 
a -particles.  ^He  reactions  not  only  add  significantly  to  the  total  positron  yield  in  flares,  but  can  also  yield  a 
positron  depth  distribution  that  peaks  higher  in  the  solar  atmosphere.  We  discuss  the  impact  these  reactions 
have  on  the  analysis  of  the  annihilation  line  observed  with  RHESSI  from  the  2002  July  23  flare.  A  significant 
contribution  from  ^He  reactions  expands  the  utility  of  the  annihilation  line  as  a  sensitive  tool  for  investigating 
the  structure  of  the  flaring  solar  atmosphere. 

Subject  headings:  acceleration  of  particles  —  nuclear  reactions,  nucleosynthesis,  abundances  —  Sun:  flares  — 
Sun:  X-rays,  gamma  rays 


1.  INTRODUCTION 

In  the  solar  photosphere,  the  abundance  of  ^He  is  believed 
to  be  quite  small.  Although  direct  spectroscopic  measurement 
of  its  photospheric  abundance  is  not  possible,  estimates  of  the 
^He/H  ratio  derived  from  solar  wind,  corona,  and  meteoritic 
measurements  range  from  1.2x10^^  (Anders  &  Grevesse 
1989)  to  3.4  X  10  ^  (Geiss  1982).  However,  in  impulsive  solar 
energetic  particle  (SEP)  events  observed  in  interplanetary  space 
(where  the  abundances  of  accelerated  particles  heavier  than 
oxygen  can  be  enhanced  by  factors  of  3-10),  ^He  can  be  en¬ 
hanced  by  several  orders  of  magnitude,  with  ^He/a  values  of 
unity  or  even  greater  (Reames,  Meyer,  &  von  Rosenvinge 
1994).  While  large  SEP  events  are  thought  to  be  accelerated  by 
shocks  in  the  solar  corona  and  interplanetary  space  (Reames 
1999,  2000),  impulsive  SEP  events  are  believed  to  originate 
low  in  the  solar  atmosphere.  Aschwanden  et  al.  (1996)  places 
the  height  of  the  impulsive  event  acceleration  region  at 
(18-20)  X  10^  km. 

Analyses  of  gamma-ray  line  flares  (e.g.,  Murphy,  Dermer,  & 
Ramaty  1987;  Share  &  Murphy  1999),  which  are  also  believed 
to  occur  in  the  lower  solar  atmosphere,  have  shown  that  the 
abundances  of  the  heavy  accelerated  particles  responsible  for 
broad  gamma-ray  lines  are  also  enhanced,  possibly  in  a  manner 
similar  to  the  impulsive  SEP  event  enhancements,  suggesting 
that  these  two  populations  may  be  accelerated  by  the  same 
process.  Observations  of  gamma-ray  lines  produced  predomi¬ 
nantly  by  accelerated  ^He,  such  as  the  0.937  MeV  line  of 
from  the  reaction  '®0(^He,/i)^^F*,  can  directly  measure  the  ^He 
abundance.  Studies  of  this  line  measured  by  the  moderate- 
resolution  Solar  Maximum  Mission  (SMM)  Gamma-Ray  Spec¬ 
trometer  (GRS)  detectors  for  several  solar  flares  (Share  & 
Murphy  1998;  Mandzhavidze,  Ramaty,  &  Kozlovsky  1997, 
1 999)  have  suggested  that  the  ^He/a  ratio  could  be  ~0. 1  or  even 
larger.  If  accelerated  ^He  is  significantly  enhanced  in  gamma-ray 


line  flares  as  these  studies  suggest,  calculations  of  yields  of  ac¬ 
celerated-particle  reactions  must  take  this  into  account. 

Recent  high  spectral  resolution  measurements  of  the 
0.511  MeV  positron  annihilation  line  (Share  et  al.  2003)  by 
RHESSI  have  focused  attention  on  the  processes  responsible 
for  the  formation  of  this  line.  Kozlovsky,  Lingenfelter,  & 
Ramaty  (1987)  provided  a  detailed  treatment  of  positron 
production  from  the  decay  of  radioactive  nuclei  produced  in 
nuclear  interactions  of  accelerated  protons  and  a -particles 
with  the  most  abundant  ambient  elements.  A  computer  code 
to  calculate  the  yields  for  various  accelerated  particle  spectra 
and  interaction  models  was  constructed  using  the  cross 
sections  given  in  that  paper.  The  code  has  been  used  in 
many  investigations  (e.g.,  Murphy  et  al.  1987)  concerning 
the  observation  of  the  0.511  MeV  annihilation  line  produced 
when  the  positrons  annihilate  with  ambient  electrons. 

In  this  paper,  we  extend  the  calculations  of  Kozlovsky  et  al. 
(1987)  to  include  the  most  important  positron-emitter- 
producing  reactions  of  accelerated  ^He.  These  ^He  reactions 
are  important  because  of  their  very  low  threshold  energies  and 
large  cross  sections  at  low  energies.  As  a  result  of  the  low 
threshold  energies,  these  reactions  may  dominate  positron- 
emitter  production  for  steep  accelerated-particle  spectra  when 
the  accelerated  particles  are  enriched  with  ^He.  As  we  show, 
the  reactions  not  only  add  significantly  to  the  total  positron- 
production  yields  but  can  also  aflhct  the  depth  distribution  of 
positron-emitter  production  in  the  solar  atmosphere.  The  in¬ 
clusion  of  ^He  reactions  is  critical  to  ensure  the  accuracy  of  the 
positron  annihilation  line  as  a  sensitive  tool  for  investigating 
both  the  flare  acceleration  process  and  the  structure  of  the 
flaring  solar  atmosphere.  We  discuss  the  impact  these  reactions 
have  on  the  analysis  (Share  et  al.  2003)  of  the  annihilation  line 
observed  with  RHESSI  from  the  2002  July  23  flare,  for  which 
two  very  different  scenarios  have  been  suggested  to  explain 
the  observed  line  profile. 
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TABLE  1 

^He-lNDUCED  Reaction  Products  and  Threshold  Energies 
(MeV  nucleon"') 

Products 


Target  n  p  “a”’’ 


'^C .  >''0  (0.49)“=  '4N(0)  "0(0)*= 

'^N .  '"F  (OA?)'  “’0  (0)  '"O  (0)°  '^N  (O)' 

“■0 .  '"NeCLl?)'  '"FU)”  '^(1.94)°  '"O  (O)' 

20Ne .  ^^MgCO)^  22]^  (0)“=  ^'NCl.l?)^  '''Ne  (0)° 

24Mg .  26si  (0.024)“=  26^1(0)“=  25^1(1.21)“=  23Mg(0)‘= 

.  203(0.165)“=  20p(0)<=  29p  (1  02)“=  (0)<= 

5<'Fe .  28^1  (0)  58Co  (0)“=  57^0  (0)  55pg  (Q) 


Note. — Threshold  energies  are  given  in  parentheses. 

““  “rf”  denotes  either  d  or  p,  n. 

“a”  denotes  either  a  or  various  combinations  of  nucleons. 

“=  Radioactive  positron  emitter. 

In  §  2  we  present  the  cross  sections  for  positron-emitter 
production  by  accelerated  ^He.  In  §  3  we  calculate  positron 
yields  for  various  accelerated  particle  spectra,  and  in  §  4  we 
relate  the  results  to  the  annihilation-line  measurements  of  the 
2002  July  23  flare  observed  by  RHESSI. 

2.  CROSS  SECTIONS 

In  Table  1,  we  list  the  products  and  corresponding  threshold 
energies  (shown  in  parentheses)  of  the  most  important  ^He 
reactions  with  the  most  abundant  elements  in  the  solar  at¬ 
mosphere.  (A  value  of  0  indicates  an  exothermic  reaction.) 
Conspicuous  features  of  this  table  are  that  most  of  the  reac¬ 
tions  are  exothermic  and  that  the  endothermic  reactions  have 
very  low  threshold  energies.  The  reason  for  these  uniquely 
low  thresholds  is  the  high  mass  excess  of  ^He  (15.8  MeV). 
This  excess  is  also  responsible  for  highly  excited  compound 
systems  in  these  reactions  and  hence  for  their  large  cross 
sections.  These  characteristics  are  the  reason  that  ^He-induced 
reactions  have  marked  advantages  for  nuclear  research  when 
only  low-energy  accelerators  (less  than  5  MeV)  are  available 
(e.g.,  Bromley  &  Almqvist  1960).  Most  of  the  products  shown 
in  Table  1  are  radioactive  positron-emitting  nuclei.  (The  ra¬ 
dioactive  products  of  Fe  reactions  mostly  decay  via  electron 
capture.)  Furthermore,  most  of  these  products  have  short 
lifetimes  (see  Kozlovsky  et  al.  1987  for  a  list  of  lifetimes)  and 
hence  can  contribute  to  the  observed  0.511  MeV  positron- 
annihilation  line  in  solar  flares,  for  which  observation  times 
are  typically  tens  of  minutes.  Kozlovsky  et  al.  (1987)  found 
these  same  radioactive  nuclei  to  be  the  most  important  con¬ 
tributors  to  positron  production  by  accelerated  protons  and 
a  -particles  in  solar  flares. 

In  the  following  sections,  we  discuss  the  cross  sections  for 
production  of  radioactive  positron-emitters  by  ^He  reactions 
with  the  various  target  nuclei  of  Table  1.  In  Figures  1-6  we 
show  the  estimated  cross  sections  of  the  various  reactions  for 
the  energy  range  of  1-100  MeV  nucleon^'.  We  discuss  the 
large  body  of  laboratory  cross  section  measurements  that  we 
use,  but  these  measurements  do  not  always  cover  the  full 
energy  range.  We  use  the  following  rules,  based  on  extensive 
studies,  to  estimate  cross  sections  for  which  measurements  are 
not  available. 

1.  For  reactions  whose  exit  channel  contains  multiple  par¬ 
ticles  of  the  form  (^He,  kpln),  where  k  and  /  are  integers  and 
k,l  >2,  the  structure  of  the  cross  section  consists  of  a  steep 


rise  from  threshold  energy  to  a  first  peak,  followed  by  a 
minimum,  and  then  a  second  rise  to  a  flat  plateau.  The  first 
peak  is  due  to  emission  of  an  a -particle  (mainly  via  a  com¬ 
pound  mechanism),  and  the  plateau  is  due  to  emission  of 
progressively  more  and  more  nucleons  (e.g.,  Gadioli  & 
Hodgson  1986;  Michel  &  Galas  1983). 

2.  For  reactions  whose  exit  channel  contains  one  particle  or 
multiple  particles  of  the  form  (^He,  kpln)  but  where  kj<2 
(e.g.,  ^He,  pn),  the  cross  sections  typically  rise  steeply  from 
threshold  energy  to  a  maximum  and  then  fall  rapidly  to  a  long 
tail. 

3.  The  steep  rise  of  the  cross  section  just  above  threshold 
energy  to  the  first  maximum  is  typical  of  reactions  between 
charged  nuclei  and  represents  the  penetration  factor  of  the 
Coulomb  barrier  (Blatt  &  Weisskopf  1952).  This  behavior  can 
be  clearly  seen  in  the  numerous  measured  cross  sections 
compiled  by  Keller  et  al.  (1973). 

4.  For  a  few  reactions  there  were  no  laboratory  mea¬ 
surements  available  to  us.  For  these  reactions,  we  have 
estimated  the  cross  sections  using  the  procedure  given  by 
Keller,  Lange,  &  Munzel  (1974,  hereafter  KLM74).  Based  on 
^1800  experimentally  determined  cross  sections  (compiled  by 
Keller  et  al.  1973),  KLM74  have  devised  an  elaborate  semi- 
empirical  procedure  for  obtaining  estimates  of  unknown  cross 
sections  for  accelerated  p,  d,  ^He,  and  a  reactions  with  nuclei 
heavier  than  fluorine. 

2.1.  Reactions  o/^He  with 

The  radioactive  positron-emitting  nuclei  produced  in  ^He 
reactions  with  are  *^N,  and  **C,  as  shown  in  Table  1. 
In  Figure  1,  we  present  the  cross  sections  for  these  reactions. 
The  data  from  2  to  8  MeV  nucleon^*  are  from  Cochran  & 
Knight  (1962),  from  0.6  to  3.3  MeV  nucleon^'  from  Hahn  & 
Ricci  (1966),  and  from  1.6  to  11  MeV  nucleon^*  from  Brill 
(1965).  At  higher  energies  the  cross  section  for  is  from 
Crandall  et  al.  (1956)  and  Aslanides  et  al.  (1981).  The  labo¬ 
ratory  measurements  of  the  cross  section  for  production  of 
*^C,  which  cover  the  entire  energy  range  of  interest  here, 
clearly  demonstrate  the  validity  of  rule  1  above. 

2.2.  Reactions  o/^He  with  ^^N 

The  radioactive  positron-emitting  nuclei  produced  in  ^He 
reactions  with  '^N  are  ^®F,  *^0,  and  *^N,  as  shown  in  Table  1. 
In  Figure  2,  we  present  the  cross  sections  for  these  reactions. 


Fig.  1. — Radioactive  positron-emitter  production  cross  sections  for  ^He 
reactions  on  '^C. 
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Fig.  2. — Radioactive  positron-emitter  production  cross  sections  for  ^He 
reactions  on 

The  data  for  *^0  from  2  to  6  MeV  nucleon^^  are  from  Hahn  & 
Ricci  (1967).  At  higher  energies,  we  followed  the  antici¬ 
pated  behavior  for  this  kind  of  reaction  according  to  the  above 
rules.  For  we  use  the  data  from  Hahn  &  Ricci  (1966) 
at  ~0.6  MeV  nucleon^'  and  from  Brill  (1965)  at  3.3-10  MeV 
nucleon^  k  Here  we  again  confirm  the  rule  concerning  the  rise 
of  the  cross  section.  The  higher  energy  rise  begins  just  beyond 
the  threshold  for  the  (^He,  ^Hen)  reaction  and  is  consistent 
with  the  explanation  that  the  rise  is  due  to  the  emission  of 
several  particles  in  the  exit  channel  rather  than  a  only. 

The  cross  section  for  production  of  is  not  available,  but 
the  (^He,  «)  reaction  is  almost  always  the  weakest  of  the 
reactions  in  Table  1  and  always  much  less  than  the  (^He,  a) 
reaction.  We  therefore  estimate  that  the  cross  section  for  *®F 
will  be  much  less  than  30  mbam,  which  is  the  value  of  the  first 
maximum  for  the  (^He,  a)  reaction  (see  Fig.  2),  so  we  do  not 
include  it  here. 

2.3.  Reactions  o/^He  with 

The  radioactive  positron-emitting  nuclei  produced  in  ^He 
reactions  with  are  **Ne,  ^^F,  '^F,  and  '^O,  as  shown  in 
Table  1.  In  Figure  3,  we  present  the  cross  sections  for  these 
reactions.  The  cross  section  for  **Ne  from  1.3  to  13  MeV 
nucleon^  ^  is  based  on  the  cross  section  measurements  by 


Fig.  3. — Radioactive  positron-emitter  production  cross  sections  for  ^He 
reactions  on  ^^0. 


TatischelT  et  al.  (2003)  for  the  reaction  ^^0(^He,  «7i  887)*^Ne. 
We  have  multiplied  this  cross  section  by  1.5  to  account  for 
additional  **Ne  production  in  the  ground  state  and  in  other 
levels  that  decay  directly  to  the  ground  state.  The  data  for  '^F 
from  0.9  to  3.2  MeV  nucleon^'  are  from  Hahn  &  Ricci 
(1966);  from  4.3  to  13.3  MeV  nucleon^*  they  are  from 
Fitschen  et  al.  (1977).  This  cross  section  includes  the  pro¬ 
duction  of  ^*F  via  the  decay  of  *^Ne.  The  data  for  from 
0.9  to  3.2  MeV  nucleon^*  are  from  Hahn  &  Ricci  (1966). 
We  obtained  a  cross  section  value  at  7  MeV  nucleon^'  from 
the  evaluation  by  Brill  (1965)  and  at  8  MeV  nucleon^*  by 
integrating  the  angular-dependent  cross  sections  of  Fuchs  & 
Oeschler  1973).  We  constructed  the  complete  cross  section 
according  to  the  rules  discussed  above.  The  data  for  *^F  are 
from  Hahn  &  Ricci  (1966)  for  2.5-3  MeV  nucleon^'.  Again, 
we  consfructed  the  complete  cross  section  according  to  the 
rules  discussed  above. 

In  addition  to  the  above  reactions  from  Table  1,  we  also 
show  in  Figure  3  the  cross  section  for  the  reaction  producing 
from  the  reaction  (^He,  2q;).  The  data  from  3.6  to  10  MeV 
nucleon^*  are  from  Brill  (1965).  At  higher  energies  we  con¬ 
structed  the  cross  section  according  to  the  above  rules  with  a 
rise  to  a  plateau. 

2.4.  Reactions  o/^He  with  ^°Ne 

All  of  the  reaction  products  of  ^He  reactions  with  ^°Ne 
shown  in  Table  1  are  radioactive  positron-emitting  nuclei. 
Because  the  half-life  of  ^^Na  is  so  long  (2.6  yr),  we  do  not 
consider  it  here.  Since  the  total  cross  sections  for  the  re¬ 
maining  nuclei  are  not  available  in  the  literature,  we  estimated 
them  by  using  the  procedure  given  by  KLM74.  The  resulting 
cross  sections  were  identical  to  cross  sections  of  the  cor¬ 
responding  reactions  of  ^He  with  ^'^Mg,  when  estimated  by  the 
same  procedure.  We  therefore  use  the  measured  ^^Mg  cross 
sections  for  the  n,  p,  d,  and  a  reactions  of  Table  1  for  the 
corresponding  ^°Ne  reactions.  The  ^^Mg  cross  sections  are 
discussed  in  the  next  section. 

2.5.  Reactions  o/^He  with  ^^Mg 

The  radioactive  positron-emitting  nuclei  produced  in  ^He 
reactions  with  ^"^Mg  are  ^®Si,  ^®A1,  ^^Al,  and  ^^Mg,  as  listed  in 
Table  1.  In  Figure  4,  we  present  their  cross  sections.  The  data 


Energy  (MeV  nucleon"'') 

Fig.  4. — Radioactive  positron-emitter  production  cross  sections  for  ^He 
reactions  on  ^^Mg.  These  cross  sections  are  also  used  as  the  cross  sections  for 
the  corresponding  n,p,  d,  and  a  reactions  of  ^*^Ne. 
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Energy  (MeV  nucleon-^) 

Fig.  5. — Radioactive  positron-emitter  production  cross  sections  for  ^He 
reactions  on  ^*Si. 


for  ^®Si,  ^^AI,  and  from  3.6  to  13.3  MeV  nucleon^'  are 
from  Frantzvog  et  al.  (1982).  At  energies  less  than  3.6  MeV 
nucleon^*  and  greater  than  13.3  MeV  nucleon^*,  we  con¬ 
structed  the  complete  cross  sections  according  to  the  rules 
discussed  above.  For  ^®A1,  we  constructed  the  cross  section 
according  to  the  procedure  given  by  KLM74.  We  assume 
that  50%  of  the  total  cross  section  produces  the  radioactive 
positron-emitter  ^®™A1. 

2.6.  Reactions  o/^He  with  ^*Si 

The  radioactive  positron-emitting  nuclei  produced  in  ^He 
reactions  with  ^^Si  are  ^°P,  ^^P,  and  ^^Si.  These  reactions 
are  listed  in  Table  1,  and  in  Figure  5  we  show  their  cross 
sections.  The  data  for  ^^P  from  3.6  to  13.3  MeV  nucleon^’  are 
from  Frantzvog  et  al.  (1982).  For  ^°P  and  ^^Si,  we  constructed 
the  cross  section  according  to  the  procedure  given  by  KLM74. 
According  to  KLM74,  the  cross  section  for  producing  is 
less  than  10  mbam,  so  we  do  not  include  it  here. 

In  addition  to  the  above  reactions  from  Table  1,  we  also 
show  in  Figure  5  the  cross  section  for  the  reaction  producing 
^^Mg  from  the  reaction  (^He,  2q;).  The  data  from  3.6  to  13.3 
MeV  nucleon^*  are  from  Frantzvog  et  al.  (1982).  At  higher 
energies,  we  constructed  the  cross  section  according  to  the 
above  rules,  with  a  rise  to  a  plateau. 

2.7.  Reactions  o/^He  with  ^®Fe 

The  only  radioactive  positron-emitting  nucleus  produced  in 
^He  reactions  with  ^®Fe  shown  in  Table  1  is  ^^Co,  in  which 
only  15%  of  the  decays  emit  a  positron.  ^^Fe  and  ^^Co  decay 
only  by  electron  capture.  Since  the  cross  section  for  total  ^^Co 
production  is  less  than  50  mbam  (Hazan  &  Blann  1965),  we 
do  not  include  it  here. 

In  addition  to  the  reactions  in  Table  1,  we  consider  pro¬ 
duction  of  ^®Co  and  ^^Ni.  A  total  of  20%  of  ^®Co  and  47%  of 
^^Ni  decay  via  positron  emission.  We  show  in  Figure  6  the 
cross  sections  for  the  production  of  ^®Co  and  ^^Ni;  the  data 
from  2  to  10  MeV  nucleon^*  are  from  Hazan  &  Blann  (1965). 

3.  POSITRON  YIELDS 

We  calculate  thick-target  positron  yields  from  accelerated 
proton,  ^He  and  a -particle  reactions  according  to  the  for¬ 
malism  given  in  Ramaty  (1986).  The  ambient  and  accelerated- 
ion  abundances  used  in  the  calculations  are  summarized  in 


Table  2.  For  the  ambient  medium  we  have  assumed  coronal 
abundances  (Reames  1995)  for  C,  N,  O,  Mg,  Al,  Si,  S,  Ca,  and 
Fe  relative  to  H  and  have  taken  Ne/0  =  0.25,  He/H  =  0.10, 
and  ^He/H  =  3x10  These  are  the  same  ambient  abun¬ 
dances  used  by  Ramaty,  Mandzhavidze,  &  Kozlovsky  (1996), 
except  for  He/H,  for  which  they  assumed  the  coronal  value  of 
0.037.  We  assume  various  values  for  accelerated  ^He/a,  and 
for  the  other  accelerated  ions  we  assume  “impulsive  flare” 
abundances  defined  by  Ramaty  et  al.  (1996),  which  are  coro¬ 
nal  (Reames  1995)  for  C,  N,  Ne,  Mg,  Al,  Si,  S,  Ca,  and  Fe 
relative  to  O  but  have  Ne/0,  Mg/0,  Si/0  and  S/0  ratios  in¬ 
creased  by  a  factor  of  3,  Fe/0  increased  by  a  factor  of  10,  and 
a /O  —  50.  This  composition  reflects  average  heavy-element 
abundance  enhancements  found  in  impulsive  SEP  events  in 
space  (Reames  et  al.  1994).  We  also  assume  ajp  =  0.5,  which 
is  at  the  maximum  of  the  range  observed  in  impulsive  SEP 
events.  Analyses  of  gamma-ray  line  flares  (Share  &  Murphy 
1997;  Mandzhavidze  et  al.  1997,  1999)  have  suggested  that 
the  large  observed  flux  of  the  a-a  line  complex  at  ~0.47  MeV 
requires  either  such  a  large  value  of  accelerated  alp  or  a 
correspondingly  large  value  of  ambient  ^He/H,  which  is 
probably  less  likely.  We  assume  a  power-law  form  for  the 
energy  spectra  of  the  accelerated  ions,  normalized  so  that  the 
number  of  accelerated  protons  greater  than  10  MeV  is  1,  and 
we  assume  that  all  species  have  the  same  spectral  index.  For 
the  proton  and  a -particle  reactions,  we  use  the  cross  sections 
given  by  Kozlovsky  et  al.  (1987),  which  include  the  pion- 
producing  reactions. 

In  Figure  7,  we  show  total  thick-target  positron  yields, 
Q+,  as  functions  of  the  accelerated-particle  power  law 
spectral  index  s  for  accelerated  ^He/a  ratios  of  0,  0.1,  0.5, 
and  1.  For  harder  spectra  (^  <  ^^3),  the  ^He  contribution  is 
negligible  even  for  He/a  =  13,  with  positron  production 
being  primarily  from  accelerated  protons  and  a  particles 
(for  the  hardest  spectra,  positrons  from  the  decay  of  pions 
produced  by  accelerated  protons  are  most  important).  For 
steeper  spectra  (s  >  ~3),  the  contribution  from  ^He  can 
dominate  (at  ^  =  5.5  and  He/a  =  13,  the  increase  is  about 
an  order  of  magnitude).  We  note  that  even  the  yields  cal¬ 
culated  here  for  He/a  =  13  (when  renormalized  to  1  proton 
>30  MeV)  are  larger  than  those  calculated  by  Kozlovsky 
et  al.  (1987),  except  for  the  hardest  spectra,  where  the  pro¬ 
duction  is  primarily  from  the  decay  of  pions  from  proton 


Fig.  6. — Radioactive  positron-emitter  production  cross  sections  for  ^Fle 
reactions  on  ^'’Fe. 
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TABLE  2 

Ambient  and  Accelerated-Ion  Compositions 


Element  Ambient  Accelerated 


H . 

1.0 

1.0 

^He . 

3.0 

X 

10-5 

a 

''He . 

0.1 

0.5 

C . 

2.96 

X 

10-'' 

4.65 

X 

10-5 

N . 

7.90 

X 

10-5 

1.24 

X 

10-5 

0 . 

6.37 

X 

10-'' 

1.00 

X 

10-5 

Ne . 

1.59 

X 

10-'' 

4.55 

X 

10-5 

Mg . 

1.25 

X 

10-'* 

5.89 

X 

10-5 

A1 . 

1.00 

X 

10-5 

1.57 

X 

10-'' 

Si . 

9.68 

X 

10-5 

4.55 

X 

10-5 

S . 

2.03 

X 

10-5 

9.56 

X 

10-'' 

Ca . 

6.75 

X 

10-5 

1.06 

X 

10-'' 

Fe . 

8.54 

X 

10-5 

1.34 

X 

10-5 

“  The  calculations  are  performed  for  various  values 
of  the  accelerated  ^He  abundance. 


interactions.  This  increase  is  due  to  the  ditferent  abundances 
assumed  here. 

Once  emitted,  positrons  either  slow  down  and  annihilate 
with  ambient  electrons  to  produce  annihilation  photons  or  es¬ 
cape  from  the  region  without  significant  annihilation  during 
typical  solar-flare  observation  times.  Positrons  annihilate  either 
directly  with  ambient  electrons  (forming  two  0.511  MeV  pho¬ 
tons  in  the  center  of  mass)  or  by  first  forming  positronium. 
Positronium  is  formed  in  one  of  two  states  depending  on 
the  relative  spins  of  the  positron  and  electron:  a  singlet  state 
(25%  of  the  time)  and  a  triplet  state  (75%  of  the  time).  Anni¬ 
hilation  from  the  singlet  state  produces  two  0.511  MeV  photons 
in  the  center  of  mass,  while  annihilation  from  the  triplet  state 
produces  three  photons  with  energies  less  than  0.511  MeV. 

Relevant  quantities  associated  with  positron  annihilation  are 
(1)  the  number  of  0.511  MeV  line  photons  produced  per 
positron  {fsn),  (2)  the  fraction  of  positrons  that  annihilate  via 
positronium  (/ps),  and  (3)  the  ratio  of  the  number  of  photons  in 
the  three-photon  continuum  to  the  number  in  the  line  (37/27). 
If  the  positronium  is  undisturbed,  its  annihilation  proceeds 
according  to  the  above  percentages  so  that  =2— 1.5/ps 
and  37/27  =  2.257ps/(2— 1.5_/ps).  However,  high  density  or 
high  temperature  can  quench  the  longer-lived  triplet  state. 


2  3  4  5  6 

Spectral  index  (s) 


Eig.  7. — Total  thick-target  positron  yield,  as  a  function  of  the  accelerated- 

particle  spectral  index  for  ^He/a  =  0,  0.1,  0.5,  and  1.0. 


either  by  ionizing  the  positronium  or  by  transforming  the 
triplet  state  to  the  singlet  state  by  spin  flipping.  This  will  tend 
to  increase /s  11  and  decrease  37/27.  If  all  of  the  triplet  state  is 
quenched  (or  if /ps  =  0),  /511  will  have  its  maximum  value  of  2 
and  37/27  =  0.  If/ps  =  1  and  the  positronium  is  undisturbed, 
/sii  will  be  0.5  and  37/27  =  4.5.  Care  must  be  taken  inter¬ 
preting  measured  values  of  these  quantities,  since  f^w  can 
appear  reduced  and  37/27  modified  because  of  (1)  finite  ob¬ 
servation  times  (since  not  all  of  the  emitters  can  decay  to 
produce  positrons  and  subsequent  annihilation  photons  due  to 
the  lifetimes  of  the  various  positron  emitters),  and  (2)  escape 
of  positrons  from  the  annihilation  region. 

The  lower  curves  of  Figure  8  show  calculated  ratios  of  the 
total  0.511  MeV  positron-annihilation  line  fluence  to  the 
summed  fluences  of  the  4.44  and  6.13  MeV  and 
deexcitation  lines,  as  functions  of  the  accelerated  particle 
spectra  index  s  for  ^He/a  =  0,  0.1,  0.5,  and  1.0.  We  assume  a 
positron  to  annihilation-line  photon  conversion  factor  fsn  = 
1.0.  If  all  positrons  annihilate  within  the  flare  observation  time 
and  there  is  no  quenching  of  the  positronium  triplet  state,  this 
factor  would  correspond  to  a  positronium  fraction  of  0.67  and 
a  37/27  ratio  of  1.5.  The  upper  curve  of  Figure  8  shows  the 
fraction  of  the  total  number  of  annihilation  photons  observable 
within  the  first  1000  s  for  instantaneous  emitter  production  at 
t  =  0,  as  a  function  of  the  accelerated  particle  spectral  index  s. 
For  very  hard  spectra  {s  <  3),  positrons  are  produced  pri¬ 
marily  via  the  decay  of  pions,  which  have  a  very  short  lifetime 
(<10^^  s),  and  nearly  all  of  the  photons  are  observable.  As  the 
spectrum  softens,  positron  production  shifts  to  radioactive 
emitters  with  longer  lifetimes,  and  fewer  photons  are  ob¬ 
servable.  For  indices  larger  than  5  there  is  a  small  dependence 
of  this  fraction  on  ^He/a,  but  the  deviation  from  the  shown 
curve  is  less  than  7%  at  5  =  6. 

In  the  next  section  we  compare  our  calculations  of  annihi¬ 
lation  line  yields  with  measured  fluences  obtained  with 
RHESSI  from  the  2002  July  23  flare.  We  also  explore  the 


Fig.  8. — Lower  curves:  Calculated  ratio  of  the  total  0.511  MeV  positron- 
annihilation  line  fluence  to  the  summed  fluences  of  the  4.44  and  6.13  MeV 
and  deexcitation  lines,  as  a  function  of  accelerated-particle  spectra  index  s 
for  ^He/a  =  0,  0.1,  0.5,  and  1.0.  The  positron  to  annihilation-line  photon 
conversion  factor  =  1.0.  Upper  curve:  Fraction  of  the  total  positron  an¬ 
nihilation  photons  observable  within  the  first  1000  s,  as  a  function  of  the 
accelerated-particle  spectral  index. 
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Fig.  9. — Calculated  ratio  of  the  0.511  MeV  positron-annihilation  line  flu- 
ence  to  the  summed  fluences  of  the  4.44  and  6.13  MeV  and  **’0  de¬ 
excitation  lines  as  a  function  of  the  accelerated-particle  spectra  index  i  for 
^He/a  =0,  0.15,  0.23,  and  1.0.  The  assumed  positron  to  annihilation-line 
photon  conversion  factor  An  is  1.0,  and  all  curves  have  been  reduced  to 
account  for  the  fraction  of  the  total  annihilation-line  fluence  observable  during 
the  23  July  flare  observation  interval.  The  ^He/avalues  of  0.15  and  0.23 
are  the  90%  and  99%  confidence  upper  limits,  respectively,  determined  from 
the  measured  ratio  of  the  0.937  MeV  **F  line  fluence  upper  limit  and  the 
~0.45  MeV  a-a  complex  fluence.  The  horizontal  band  indicates  the  ratio 
(±1  a)  observed  by  RHESSl  from  the  2002  July  23  flare.  The  vertical  dotted 
lines  indicate  the  allowed  range  (±1  a)  of  spectral  index  derived  from  the 
observed  fluence  ratio  of  the  2.223  MeV  neutron-capture  line  to  the  summed 
4.44  and  6.13  MeV  and  *'’0  deexcitation  lines. 


consequences  of  these  calculations  for  the  high  spectral  res¬ 
olution  measurements  of  the  annihilation  line  shape. 

4.  APPLICATION  TO  THE  2002  JULY  23  SOLAR  FLARE 

Figure  9  shows  calculated  ratios  of  the  0.511  MeV  annihi¬ 
lation-line  fluence  to  the  summed  4.44  MeV  and  6.13  MeV 
*®0  deexcitation-line  fluences,  as  functions  of  the  accelerated- 
particle  spectral  index  s  for  the  extreme  cases  of  ^He/ct  =  0 
and  1.0.  We  have  taken  into  account  the  fraction  of  the 
annihilation-line  fluence  expected  to  be  observed  during  the 
960  s  observation  time  of  the  July  23  flare,  assuming  that 
the  positron-emitter-production  time  history  is  identical  to  the 
nuclear  line-emission  time  history  for  this  flare  (Share  et  al. 
2003).  Upper  limits  for  the  accelerated  ^He/a  abundance  ratio 
for  this  flare,  obtained  from  a  comparison  of  the  mea¬ 
sured  upper  limits  for  the  ^He-produced  0.937  MeV  line  flu¬ 
ence  and  the  fluence  of  the  a -a  complex,  are  0.15  with  90% 
confidence  and  0.23  with  99%  confidence.  These  values  are 
comparable  to  the  SMM  19-flare  average  for  ^He/a  of  0.1 
(Share  &  Murphy  1998).  Calculated  fluence  ratios  for  these 
values  of  ^He/a  are  also  shown  in  Figure  9.  For  all  of  the 
curves  of  Figure  9, /sn  was  assumed  to  be  1,  which  corre¬ 
sponds  to  a  positronium  fraction of  0.67  and  a  37/27  ratio 
of  1.5.  This  value  for  37/27  is  consistent  with  the  measured 
99%  confidence  upper  limit  of  3.3  (Share  et  al.  2003).  We  note 
that  there  is  some  dependency  of  the  calculated  ratio  on  the 
assumed  accelerated  alp  abundance  ratio.  An  estimate  of  this 
abundance  ratio  for  the  23  July  flare  can  be  obtained  by 
comparing  calculated  fluence  ratios  of  the  a-a  complex  and 
the  summed  4.44  and  6.13  MeV  lines  with  the  measured  ratio 
of  2.2  ±  0.9  (Share  et  al.  2003;  Smith  et  al.  2003).  The 
measured  ratio  provides  a  1  ct  lower  limit  of  ~0.35,  which  is 
consistent  with  our  assumed  value  of  0.5. 


Also  shown  in  Figure  9  is  the  measured  line  ratio  for  this 
flare,  1.35  ±  0.5  (Share  et  al.  2003;  Smith  et  al.  2003),  indi¬ 
cated  by  dashed  lines.  The  vertical  dotted  lines  indicate  the 
range  of  spectral  index  deduced  (2004,  in  preparation)  from 
the  measured  ratio  of  the  2.223  MeV  neutron-capture  line 
fluence  to  the  summed  fluences  of  the  4.44  MeV  and 
6.13  MeV  lines.  The  calculated  line  ratio  is  consis¬ 
tent  with  the  measured  ratio  for  ^He/a  ratios  up  to  at  least  the 
99%  confidence  value  of  0.23.  We  note  that  if  the  alp  ratio 
were  as  low  as  0.1,  the  curves  of  Figures  8  and  10  in  the 
spectral  range  of  interest  would  be  reduced  to  about  70%  of 
their  values  for  a/p  =  0.5,  and  the  curves  for  ^He/a  ratios 
up  to  0.23  would  still  be  consistent  with  the  measured  line- 
fluence  ratio.  We  also  note  that  even  if  as  many  as  50%  of  the 
positrons  were  to  escape  from  the  region  without  apprecia¬ 
ble  annihilation  during  the  observation  time,  the  calculated 
ratio  for  ^He/a  =  0.23  would  still  be  consistent  with  the 
measurement. 

The  RHESSl  high-resolution  Ge  detectors  resolved  the 
0.511  MeV  positron-annihilation  line  observed  in  the  2002 
July  23  solar  flare  (Share  et  al.  2003).  When  fit  with  a 
Gaussian  line  shape,  the  line  was  found  to  be  broad  with  a 
full  width  at  half  maximum  (FWHM)  of  8.1  ±  1.1  keV. 
Share  et  al.  (2003)  found  that  the  line  shape  was  consistent 
with  two  very  different  production  scenarios.  The  line  could 
have  been  formed  in  a  quiet  solar  atmosphere  at  a  temper¬ 
ature  of  ~6000  K.  Under  these  conditions,  the  line  is  com¬ 
posed  of  both  a  narrow  and  a  broad  component.  The  narrow 
component  (~1.5  keV)  is  produced  by  annihilation  of  ther- 
malized  positrons  with  bound  electrons,  while  the  broad  com¬ 
ponent  (~7  keV)  results  from  positronium  formation  by 
charge  exchange  in  flight  with  hydrogen.  For  a  narrow  range 
of  temperatures,  the  broad  component  is  sufficiently  strong  to 
produce  an  effective  “broad”  line.  Alternatively,  the  line 
could  be  a  thermally  broadened  Gaussian,  whose  width  cor¬ 
responds  to  a  temperature  of  (4-7)  x  10^  K. 

Both  of  these  scenarios  present  difficulties.  The  difficulty 
with  the  first  scenario  is  that  it  requires  an  exceptionally  fine 


H  density  (cm-^) 

10"  10'^  10'^  10"  10'^  10'®  10'^  10'® 


Depth  d  (km) 

Fig.  10. — Fraction  of  production  occurring  deeper  than  a  given  depth  for 
two  reactions  whose  cross  sections  have  a  low-energy  threshold  (1  MeV 
nucleon"’,  solid  line)  and  a  high-energy  threshold  (25  MeV  nucleon"’,  dotted 
line).  The  dashed  lines  indicate  the  50%  points  on  each  curve.  The  assumed 
atmospheric  model  is  Avrett  (1981).  The  bottom  axis  is  depth  in  the  atmo¬ 
sphere  (km)  and  the  top  axis  is  the  corresponding  hydrogen  number  density 
(cm"^). 
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tuning  of  the  annihilation  region,  since  only  for  a  very  narrow 
range  of  temperatures  (approximately  5650-6270  K)  is  the 
broad  component  sufficiently  strong  to  produce  an  acceptable 
fit  to  the  line  shape.  At  temperatures  immediately  above  and 
below  this  range,  the  line  becomes  very  narrow.  For  the  sec¬ 
ond  scenario,  calculations  (Hua,  Ramaty,  &  Lingenfelter  1989) 
for  a  relatively  hard  accelerated-particle  spectrum  have  shown 
that  the  interaction  site  of  nuclear  reactions  similar  to  those 
responsible  for  positron-emitter  production  by  protons  and 
a -particles  is  expected  to  be  at  hydrogen  densities  of  10^^ 
to  10'^  cm^^.  Current  models  of  the  solar  atmosphere  indi¬ 
cate  that  the  temperature  at  such  high  densities  is  much  less 
than  (4-7)  X  10^  K. 

The  new  reactions  of  positron-emitter  production  by  ac¬ 
celerated  ^He  presented  here  may  help  make  the  second 
scenario  more  plausible.  The  solar  atmosphere  during  flares  is 
known  to  be  significantly  altered  from  quiet-Sun  conditions 
by  processes  such  as  heating  and  mass  motion.  There  is  ev¬ 
idence  for  high  temperatures  at  high  densities,  for  example,  in 
the  dramatic  enhancement  over  quiet-Sun  values  of  C  iv  and 
Si  IV  line  emission  in  the  transition  region,  as  noted  by  Brekke 
et  al.  (1996).  As  another  example,  Raymond  et  al.  (2003) 
observed  large  emission  measures  of  10^  K  gas  during  the 
2002  July  23  flare.  Although  temperatures  as  high  as  (4-7)  x 
10^  K  at  the  high  densities  10*^-10*^  cm^^  are  still  unlikely, 
if  a  significant  fraction  of  the  positrons  were  formed  at  sub¬ 
stantially  lower  densities,  as  will  be  shown  below,  the  high 
temperatures  implied  by  the  annihilation  line  width  would  be 
much  more  plausible  there. 

As  discussed  above,  the  ^He  reactions  are  unique  in 
that  their  cross  sections  are  both  large  and  have  very  low 
threshold  energies  compared  with  corresponding  proton-  and 
a -particle-induced  reactions.  For  a  power-law  spectrum  of 
accelerated  particles  losing  energy  in  an  ambient  medium, 
such  low-threshold  reactions  are  produced  predominantly 
by  particles  with  low  initial  energies.  These  particles  have 
shorter  ranges  than  higher  energy  particles,  so  the  reactions 
tend  to  occur  at  higher  elevations  (lower  densities)  as  the 
particles  are  moving  downward  in  the  solar  atmosphere.  We 
demonstrate  this  by  calculating  the  depth  distribution  of  nu¬ 
clear  reactions  having  a  low  threshold  energy  (~1  MeV 
nucleon^  \  similar  to  those  of  typical  positron-emitter  reac¬ 
tions  by  accelerated  ^He)  and  nuclear  reactions  having  a 
high  threshold  energy  (^25  MeV  nucleon^*,  similar  to  those 
of  typical  positron-emitter  spallation  reactions  by  accelerated 
protons). 

To  perform  these  calculations,  we  use  the  loop  model  of 
Hua  et  al.  (1989),  which  we  briefly  describe.  The  model 
includes  energy  losses  due  to  Coulomb  collisions,  removal  by 
nuclear  reactions,  magnetic  mirroring  in  the  convergent  flux 
tube,  and  MHD  pitch-angle  scattering  (PAS)  in  the  corona  (see 
Hua  et  al.  1989  for  details).  The  accelerated  ions  are  released 
isotropically  at  the  top  of  a  magnetic  loop  and  are  followed 
until  they  interact  or  thermalize,  usually  near  the  loop  foot- 
points  in  the  chromosphere  or  upper  photosphere.  The  model 
consists  of  a  semicircular  coronal  portion  and  two  straight 
portions  extending  vertically  from  the  transition  region 
through  the  chromosphere  into  the  photosphere.  Below  the 
transition  region,  the  magnetic-field  strength  is  assumed  pro¬ 
portional  to  a  power  6  of  the  pressure  (Zweibel  &  Haber 
1983).  PAS  can  be  characterized  by  its  mean  free  path  for 
isotropization.  A,  here  expressed  by  A,  the  ratio  of  A  to  the 
loop  half-length.  The  level  of  PAS  directly  affects  the  angular 
distribution  of  the  accelerated  particles  when  they  interact 


with  the  solar  atmosphere.  Using  gamma-ray  line  data  from 
several  SMM  flares.  Share  et  al.  (2002)  showed  that  the 
measured  deexcitation-line  Doppler  shifts  imply  interacting- 
ion  angular  distributions  that  are  inconsistent  with  no  PAS 
(a  oo),  obtaining  better  fits  with  values  of  A  =  300  or  less. 
Here  we  assume  a  level  of  PAS  given  by  A  =  300.  The  at¬ 
mospheric  model  is  the  sunspot  model  of  Avrett  (1981),  and 
the  accelerated-particle  power-law  spectral  index  was  assumed 
to  be  4,  typical  of  gamma-ray  line  flares  (Ramaty  et  al.  1996). 
We  assume  a  loop  length  of  1.15x10®  cm  and  6  =  0.2, 
values  found  by  Hua  et  al.  (1989)  to  provide  good  fits  to 
deexcitation-line  decay  times  in  the  1980  June  21  flare. 

For  the  two  threshold  energies,  1  and  25  MeV  nucleon^*, 
we  show  in  Figure  10  the  fraction  of  reactions  occurring 
deeper  than  a  given  depth  in  the  solar  atmosphere,  with  dashed 
lines  indicating  the  50%  points.  Reactions  with  a  low 
threshold  occur  at  significantly  higher  elevations,  with  a  dif¬ 
ference  in  density  of  more  than  two  orders  of  magnitude  for 
this  model  atmosphere.  Thus,  positron-emitter  production 
from  the  low-threshold  ^He  reactions  would  occur  higher  in 
the  chromosphere  and  at  lower  densities  than  production  from 
proton  and  a -particle  reactions. 

If  the  temperature  of  this  lower-density  region  were  suffi¬ 
ciently  high  during  the  flare,  these  positrons  would  add  a  broad 
thermal  component  to  the  line  profile.  However,  Figure  9 
shows  that  for  the  July  23  ^He/a  upper  limit  of  0.23,  the  ^He 
reactions  would  only  contribute  ^30%  to  the  0.511  MeV  line. 
Even  if  these  positrons  annihilated  in  a  high-temperature  re¬ 
gion,  their  contribution  to  the  line  would  not  be  sufficient  to 
produce  a  line  shape  consistent  with  the  width  measurement 
for  this  flare. 

5.  SUMMARY 

We  have  shown  that  accelerated  ^He  can  produce  a  signif¬ 
icant  number  of  positrons  in  solar  flares  with  steep  acceler¬ 
ated-particle  spectra,  particularly  if  the  ^He/a  ratio  is  >1,  as 
observed  in  some  impulsive  SEP  events  in  interplanetary 
space.  Because  of  their  low  threshold  energies,  these  ^He 
reactions  tend  to  occur  at  lower  densities  in  the  solar  atmo¬ 
sphere  than  the  reactions  of  protons  and  a -particles.  If  the 
temperature  at  these  lower  densities  is  higher  than  inferred 
from  current  atmospheric  models,  such  positrons  would  add  a 
broad  thermal  component  to  the  line  profile.  However,  for  the 
2002  July  23  flare  observed  by  RHESSl,  the  measured  99% 
confidence  upper  limit  for  the  ^He/a  ratio,  0.23,  is  probably 
not  sufficient  to  fully  explain  the  observed  line  shape  for  this 
flare. 

A  full  understanding  of  the  annihilation  line  will  require 
a  comprehensive,  self-consistent  calculation  of  annihilation¬ 
line  production.  The  calculations  must  include  the  transport  of 
the  accelerated  particles  and  depth  distribution  of  positron- 
emitter  production,  transport  of  the  emitted  positrons,  and 
all  possible  processes  by  which  positrons  annihilate  (direct 
annihilation,  positronium  formation,  quenching,  etc.).  All  of 
these  calculations  must  be  performed  within  the  context  of  a 
realistic  magnetic-loop  model  assuming  various  solar  atmos¬ 
pheres.  This  comprehensive  study  will  allow  us  to  construct 
an  annihilation  line  profile  that  is  composed  of  contributions 
from  a  range  of  physical  conditions  determined  by  the  depth 
distribution  of  positron  production,  transport  and  annihila¬ 
tion.  All  of  the  processes  involved  in  the  formation  of  the 
positron-annihilation  line  are  sensitive  to  the  physical  state 
of  the  flaring  loop  such  as  temperature,  density,  ionization, 
and  magnetic-field  configuration.  This  line,  whose  previous 
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theoretical  treatments  assumed  only  simplified  conditions, 
will  become  an  exceptionally  sensitive  tool  for  exploring  the 
physical  conditions  of  the  solar  flare  region. 
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